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The pion electromagnetic form factor is calculated within the QCD light-cone sum rule method 
and using a renormalon model for the higher twist distribution amplitudes (DAs). The theoretical 
predictions are compared with the experimental data and constraints on the pion leading and twist-4 
DAs are extracted. An upper bound on the twist-4 contribution to the form factor and estimates 
of effects due to higher conformal spins in the pion twist-4 DAs are obtained. 
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I. INTRODUCTION 



The leading and higher twist distribution amplitudes (DAs) of hadrons are important ingredients in investigation of 
various exclusive processes within QCD 0, Q| • The leading twist DAs appear in the QCD factorization formulas and 

■ describe exclusive processes with the leading power accuracy. They correspond to parton configurations of hadrons 
^ ' with a minimal number of constituents. The higher twist DAs are essential in computing different power-suppressed 

, corrections, which emerge due to parton virtuality, transverse momentum, contributions of higher Fock states with a 
nI " nonminimal number of hadron constituents. 

The traditional method for the description of DA is founded on the conformal symmetry of the QCD Lagrangian 
0. Within this approach the leading and higher twist DAs are expanded over the conformal spin. It is important 
ly-j \ that any parametrization of DA based on a truncated conformal expansion is consistent with the QCD equations of 
motion (EOM) and is preserved by the QCD evolution to the leading logarithmic accuracy In |2j. Therefore, the 
conformal expansion provides a practical framework for modeling of hadron DAs 0, 0, IS S Ifl U an d is widely used 

■ for investigation of numerous exclusive processes in QCD. 

I ' Because of the increasing number of parameters at higher conformal spins and practical difficulties in phenomeno- 
©H. logical applications, one has to restrict one's self by only the first few terms in the conformal expansion of DAs. As a 
result, the contributions of higher conformal spins to DAs in the existing calculations are neglected. At the same time, 
the suppression of higher spin contributions and the convergence of conformal expansion at present experimentally 
accessible energy regimes is by no means obvious and may be wrong. Therefore, one needs to draw new approaches 
to clarify this problem. 

The renormalon model proposed recently in Refs. pursues to test precisely this issue, that is to set a plausible 

upper bound for the possible contributions of higher conformal spins that so far escaped attention. The renormalon 
approach employs the assumption that the infrared (IR) renormalons in the leading twist coefficient functions should 
cancel the ultraviolet (UV) renormalons in the matrix elements of twist-4 operators in a relevant operator product 
expansion. Such cancellation was proved by explicit calculations in the case of the simple exclusive amplitude involving 
pseudoscalar and vector mesons |Tl| . It turned out that this is enough to obtain the full set of two- and three-particle 
twist-4 DAs of pseudoscalar and vector mesons in terms of the leading twist DAs. It is remarkable that the set of 
twist-4 DAs depend only on one new parameter, which can be related to the matrix element of some local operator 
(see Sec. II ) and estimated using the QCD sum rule. In other words, the twist-2 and twist-4 DAs of pseudoscalar 
and vector mesons can be determined using the same set of parameters that considerably restricts a freedom in the 
choice of DAs, increasing, at the same time, the predictive power and reliability of QCD results. 

The renormalon calculus was employed in Ref. 11], where the pion and p-meson twist-4 DAs were constructed. In 
the calculations the mesons asymptotic DAs were used. A generic feature of the renormalon model is that it predicts 
higher twist distributions that are larger at the end points compared to the asymptotic distributions, and are expected 
to give rise to larger higher twist effects in exclusive reactions. The main purpose of this work is to test this idea on 
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example of the pion electromagnetic form factor (FF), that is to set an upper bound on possible twist-4 contributions 
to FF. To this end, we extend results of Ref. and compute the pion higher twist DAs using the leading twist DA 
with two nonasymptotic terms. We apply our predictions for studying the pion form factor within the QCD light-cone 
sum rule (LCSR) method and extract constraints on the input parameters 62(^0) an d ^4(^0) a ^ the normalization 
scale Mo = 1 GeV 2 . 

This paper is structured as follows: In Sec. II we define the two- and three-particle twist-4 DAs of the pion and 
calculate them within the renormalon approach. In Sec. Ill general expressions for the FF F n (Q 2 ) in the context of 
the QCD LCSR method with twist-6 accuracy are presented. In Sec. IV we confront our predictions with the available 
data on Fk(Q 2 ) and by this way model the pion DAs. Section V is reserved for concluding remarks. 

II. HIGHER TWIST DAS OF THE PION 

The light-cone two-particle distribution amplitudes of the pion are defined through the light-cone expansion of the 
matrix element, 

(0 \d(x 2 )~fv~f5 [x 2 ,Xi]u(xi)\ TT + (j))) = 



= if w p v / due~ lupxi - zupX2 
Jo 



+ iU(A v (pA)-p v A 2 ) I due- iu ^- iupx * ^(u, mf) + 0(A 4 ) 



(2.1) 



where (p^(u, n%) — f^u, is the leading twist DA of the pion, whereas tp^(u, /J-%), ip 2 4 \u, Mf) are two-particle 
twist-4 DAs. We use the notation [x 2 , x\] for the Wilson line connecting the points x\ and x 2 , 



[x2,Xx] = Pexp 



-ig [ dtA li A' i (x2 + tA) 
Jo 



(2.2) 



In Eqs. (|2.1|) and l|2.2|) A = x\ - x 2 and u = 1 — u. 

Apart from the two-particle DAs there exist the three-particle twist-4 DAs involving an extra gluon field, which we 
define in the form 0, [Tl| 

(0 \d(-z) [-z,vz]j u j 5 gG fip (vz) [vz, z]u{z)\tt + (p)) 



U [ D ai e- wz(ai - a2+a3v) { — 
J [pz 



(PpZ p -p p z^»(ai,a2,a 3 ) 



where the longitudinal momentum fraction of the gluon is and the integration measure is defined as 



(2.3) 



J Dai = J dctida 2 da?,5(l — a% — a 2 — a 3 ). (2-4) 

The other pair of DAs that can be obtained from Eq. I|2.3|l after the replacement 75 G^p —> iG pp — k^ pa ^Gai3: 

'0 d(-z) [-z,vz]j v igGpp(vz) [vz, z] u(z) ir + {p)] 



/ Date 



-ipz(a\ — v) J Pv 

pz 



{PpZp - p p z fl )m\\(ai,a>2,a 3 ) 
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pz / \ pz 
One more three-particle DA S„-(ai) is introduced through the formula 

(0 \d(-z) [~z,vz]j^ 5 gD a G ap (vz) [vz, z] u(z)\ n + (p)) 



(2.5) 



(2.6) 



In this work we do not consider twist-4 four quark operators and corresponding DAs. 

Because of the EOM the two-particle DAs ip^\u, fi F ), ip^\u, fi F ) can be expressed in terms of the three-particle 
ones. Namely, from exact operator identities 12], which relate integrals over v of the quark-gluon-antiquark operator 
in Eq. I|2.3|) to derivatives of the quark-antiquark operator (I2.1|l , it follows that 

ip 2 4 \u)= dv dai / da 2 — [2$j_ - $11] a 2 , a 3 ), 

JO JQ JQ a 3 



da 



UCti — UCt2 
2 

ai 



[2$ ± - $||] (ai,a 2 ,a 3 ), 



(2.7) 



where = 1 — a.\ — 

The standard method to handle meson DAs is modeling them employing the conformal expansion. Then for the 
leading twist pion DA we get 



(f w (u, flp) = ipasy(u) 



1+ £ M/4)c£' a (t 



(2.8) 



Here ip asy (u) is the pion asymptotic DA 



¥>asy(u) = 6uu, 



and CfJ 2 {^) are the Gegenbauer polynomials. The functions b n ([i F ) determine the evolution of ^{u, /J/p) on the 
factorization scale /j, F , 



&n(/4") = M/^o) 



as(^ 



.«s(Mo) 



■y n //3 



n+1 



(n+l)(ra + 2) 



3=2 J 



(2.9) 



In the above, 7„ are the anomalous dimensions and [Iq is the normalization scale. The expansion over the conformal 
spin j can also be performed for the higher twist DAs |7J, |8|, lid, ll^ • 

The renormalon approach to the higher twist DAs is based on another idea. To explain principle points of the 
renormalon approach and derive relations between the pion twist two and four DAs in Ref . 1 1] , the authors considered 
the gauge-invariant time-ordered product of two quark currents at a small light-cone separation, 



(0 \T{d(x2)"fv"f5 [3:2,0:1] 7T + (p)) 



= iU f due~ iupxi -" ip] 
Jo 

with j A 2 j <C 1/A 2 and A 2 < playing the role of the hard scale and jj? being the ultraviolet renormalization scale. 
This martix element is parametrized in terms of two Lorentz-invariant amplitudes Gi(u, A 2 , /1 2 ) and G 2 (u, A 2 , /i 2 ), 
which in the light-cone limit A 2 — > with pA fixed have the expansions 



fnK/.A-. //-)/>,. + (7.i„.A-. /r) ( ^A„-p, 



(2.10) 



Gi(u, A 2 ) = Cf > ® ^ + A 2 £ ® A 4) + 0(A 4 ). 



(2.11) 
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Here cf ) are the coefficient functions and tpw are the pion DAs, t refers to twist, and the summation runs over all 
contributions for a given twist. 

Considering in l|2.11[l the leading twist coefRcient functions to all orders of as(fi 2 ) and the twist-4 contribution to 
the leading order one gets 

Gi(m,A 2 ) = [I + cia s + c 2 a 2 s + ....] <g> <^> (2) + AV[ 4) («, 



G 2 {u,A 2 ) = [c ias + c 2 a 2 s + ....] ®tpW + A 2 cp¥\u,fi? F ) 



(2.12) 



Calculation of the leading twist coefficient functions to all orders using the running coupling method gives rise to IR 
renormalon ambiguities in the amplitudes Gi(u, A 2 ). These ambiguities are expressable in terms of the pion leading 
twist DA. The twist-4 DAs (|2.3|) . (|2.5() . and l|2.6|) contain UV renormalon divergences, which were employed in Ref. 

to compute UV renormalon ambiguities in the pion two-particle twist-4 DAs <pf (u, [i 2 F ). These UV renormalon 
ambiguities cancel the IR renormalon ones in the sum of the different twists ()2.12(l . in the same way as the logarithmic 
scale dependence is cancelled between matrix elements and coefficient functions for a given twist. As a result, the 
structure functions Gi(u, A 2 ) are unambiguous to the twist-4 accuracy. The idea of the renormalon model for the 
pion twist-4 DAs is to define them by taking the functional form of the corresponding UV renormalon ambiguities 
and replacing the overall normalization constant by a suitable nonperturbative parameter. By this way one obtains 
the following relations between the DAs of the pion: 



s;2 r 



$_i_(ai, a 2 , 0:3) 



1 — Ul 1 — OL 2 



$||(ai,a2,a 3 ) 



Q2^tt(Qi) _ ai^jr(a 2 ) 
(l-a x ) 2 (l-a 2 ) 2 



1 — a.\ 1 — a 2 



*ll(ai,a 2 ,a 3 ) 



5 1 
3 



-7r(ai, a 2 , a 3 ) 



2f5 2 ["02^(01) _ a\<£-K (a 2 ) 
3 1 — ol\ 1 — a 2 



(2.13) 



As is seen, the renormalon model for the set of twist-4 DAs depends only on one free parameter S 2 . It is related to 
the matrix element of the local operator 



5 2 {nl) ~ 0.2 GeV 2 



(2.14) 



and estimated from various 2-point QCD sum rules 0] . 

In the case of the asymptotic DA, the pion twist-4 DAs were computed in Ref. ^lj. In this paper we apply the 
results of Ref. ;llj to a more general situation. To this end, we rewrite the leading twist DA l|2.8|l in the form 



lfi n (u, /if,) = ipasy(u) 2_] K n (p 2 p)u n 



n=0 



(2.15) 



This form is more suitable for calculations and leads to compact expressions for the higher twist DAs. The DA 
(^(u, n 2 F ) can also be expanded over u with the same coefficients K n {^j? F ) and, hence, Eq. I|2.15|) preserves the 
symmetry of the distribution amplitude under the replacement u <-> u, even if this is not explicitly seen from Eq. 
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For DAs containing two nonasymptotic terms C^ 2 (u — u) and C^ 2 {u — u), the sum i|2.15[l runs over n — 0, 1, ..4 
and the coefficients K n {^i 2 F ) are given by the equalities 

Xo(/x|) = 1 + 66 2 (/4) + 156 4 (/4), ^i(Mf) = -30 [& 2 (/4) + 7& 4 (/&)] » 
K 2 ( M |) = 30 [6 2 (/x|) + 286 4 (/4)] . ^04) - -12606 4 (4), 



X 4 (/4) = 6306 4 (/4)- 



(2.16) 



Calculation of the three-particle DAs (|2.13|) is straightforward. The two-particle DAs ipf^(u, jxV) and ip^ (it, fj, F ) have 
the form 

4 



n=0 



n=0 



Their components </?^(it) and ^> 2 (m) are given by the following expressions: 

( Po( u ) — 3 2 |^ [ U1 " — Li 2 (u)] + u [In it — Li 2 (u)] — uu + — 1 



(2.17) 



It U' . 

l + --y ]ln¥-Li 2 (u) 



it it' . 
1+ 2 -y )ln«-Li 2 (u) 



5 1 , , 71 

6 2 6 



</3 2 (u) = 5 2 I u 



1 + u — -u 2 ) In it — Li 2 (u) 



1 + u — —u 2 ) lnu — Li 2 (u) 



2 _ 5 2—2 71 

-— H — u u H — , . 

3 4 6 ' 



2 J u 



. 3 7 o 1 ^ 1 4 \ . 

1 + -U--U +-u - —it )lnM-Li 2 (u) 



3 7o ^ ^ 4 \ r * 

1 + -w — -it + —u — —u mit — Li 2 (u) 
2 6 4 10 



31 257 n_ 2 1 , , 7T 2 

— m H u u u it H f . 

60 120 3 6 ' 



<p\(u) = S 2 



2 J T7 



113 3 

1 + 2u - — u 2 + -u 3 - — u 4 I lnu ~ Li 2 (u) 
6 4 10 
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FIG. 1: The components of the two-particle twist-4 DAs <Pn( u ) ( a ) an d — L Pn{ u ) (b) as functions of u. The normalization 
constant is chosen equal to 5 2 = 1. 
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and 



23 47 9 9 61 o o 7r 
-uu H u u u u + 



60 15 45 



6 I ' 



</3q(u) = 5 2 [it 2 In M + u 2 In u + mm] , 



(2.18) 



¥>?(«) = <* 2 

vl(u) = S 2 
<p 2 (u) = 6 2 



u 2 In u + ii 2 In U + HU+ 



u 2 In u + m 2 In u + uu H — m 2 u 2 
6 



2 lntt + u 2 lnw + uu + y^V — -u 3 u ' 



u 2 In it 



u 2 \nu + uu - 



77 2 _2 13 3_ 3 
60 UM "30 UU 



(2.19) 



where Li a (x) = Y^Li x "/ na - 

The functions ¥>„(it) and are shown in Fig.^ As is seen, the shapes of the functions <p 2 (w) are identical to 

each other, difference being only in their normalization. On the contrary, the functions tp\(u) and f\(u) differ from 
Vo(i 2)( u ) a ^ so m their shapes and have minima at the point uq = 1/2 . With the constant 6 2 (pq) being fixed from 

the QCD sum rule, the twist-4 DAs ip^ (u, ^jl 2 f ) and ip^iu, as well as ones that are given by Eq. I|2.13|l depend 
only on the parameters &2(Mo), ^(Mo)- I n other words, in the framework of the renormalon approach the twist-4 DAs 
of the pion are determined by the function (^(it, /u^) unambiguously. 



III. THE PION ELECTROMAGNETIC FORM FACTOR WITHIN THE QCD LCSR METHOD 

In this section we apply the twist-4 DAs for calculation of the pion electromagnetic form factor. We use the QCD 
LCSR method, which is one of the powerful tools to estimate nonperturbative components of exclusive quanities 16]. 
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The LCSR expression for the pion electromagnetic FF was derived in Refs. |l7lll8|. It was reanalyzed recently in Ref. 
|19|. where a sign error in the previous calculation of the twist-4 contribution to FF was corrected. Our approach 
to the twist-4 term leads to further improvement of the prediction for FF, because the twist-4 DAs obtained in the 
previous section encompass contributions arising from higher conformal spins. 

It is worth noting that the renormalon technique was successfully employed for studying the light mesons elec- 
tromagnetic and transition FFs [2£j, El H2 • I n ^ ne works [2£j, El , the power-suppressed corrections to these FFs 
were found using the running coupling method. The latter leads to Borel resummed hard-scattering amplitudes of 
the relevant subprocesses and necessitate calculation of the QCD factorization formulas applying the principal value 
prescription. In Refs. [20L l2l| it was demonstrated that the running coupling method allows one to take into account 
both the hard and soft components of the FFs. 

The LCSR method is based on the analysis of the correlation function 

T flv (p,q)=i J d 4 xe^ x (0\T {jl(0)f u m (x)}\n+(p)) , (3.1) 

where j'fj = dj^^u and H m = e u u'ji,u + eddj^d is the quark electromagnetic current. The contribution of the pion 
intermediate state is given by 

T^ip, q) = 2iU(p - qj^F^Q 2 )— \ -. (3.2) 

Here /„• is the pion decay constant and F 7T (Q 2 ) is the pion electromagnetic FF. Because p 2 = m 2 and q 2 = —Q 2 , the 
correlation function (|3.1I) actually depends on one invariant variable s — (p — q) 2 . For large negative values of s and 
q 2 , this correlator can be computed in QCD. In the QCD sum rule method by matchig between the dispersion relation 
in terms of contributions of hadronic states and the QCD calculation at Euclidean momenta, one can estimate the 
hadronic quantities under consideration, in our case, the pion FF Ft T {Q 2 ). This is the common idea sharing by QCD 
sum rule methods, the difference being in approaches to compute the correlation function (|3.1J) within QCD. 

When the q 2 and (p — q) 2 are spacelike and large, the correlation function can be expanded near the light-cone in 
terms of the pion DA of increasing twist. As a result contributions to F 7T (Q 2 ) coming from the pion DAs of different 
twists can be found. The leading twist (twist-2) light-cone sum rule for F 1T {Q 2 ) is (hereafter m 2 — 0) ^3] 



duipn^u, fi F ) exp 



uQ 2 



uM 2 



where 



u 



Q 2 



s + Q 2 



(3.3) 



(3.4) 



In Eqs. (|3.3|) and 13. 4|) sq is the duality interval; M 2 is the Borel variable. 

The accuracy of the LCSR (|3.3() was improved by calculating O(as) correction to the twist-2 part, as well as 
including into consideration twist-4 and twist-6 contributions [l3,[l3. Finally, the F^(Q 2 ) takes the following form: 



F^Q 2 ) = F^{Q 2 ) + Fi 2 -^\Q 2 ) + F^(Q 2 ) + F^(Q 2 ). 



(3.5) 



The details of calculations and the explicit expression for Fi 2 ' (Q 2 ) can be found in Ref. Here we only remark 
that, namely, this contribution provides the standard QCD asymptotics ~ as/Q 2 of the form factor. 



The twist-4 term Fi 4 \Q 2 ) is given by expression 



F^\Q 2 )= j 

J Uy 



tp4,(u,H 2 F ) 
du rrrr— exp 



uM 2 



uQ 2 



uM 2 



u a (p 4 (u , n 2 F ) _ so/M 2 

Q 2 



where 



2u 



d (4)/ 2x - d 2 (4), 2 N 

— <p\ (u, n F ) ~ u-^ip\ (u, Hp) 



(3.6) 



(3.7) 



The difference between Eq. (|3.7f) and the relevant formula in Ref. 19] is connected with the definition of the distribution 
amplitude (p^\u, fi F ). In fact, the twist-4 DA gi^iu^p) used in Ref. [Tiil ] can be written in terms of ip^ (w, fi F ) 
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FIG. 2: The dependence of the light-cone sum rule on the Borel parameter. The asymptotic DA is used. For the solid curve 
Q 2 = 1 GeV 2 , for the dashed curve Q 2 = 4GeV 2 , and for the dot-dashed one Q 2 = 10 GeV 2 . 



The factorizable twist-6 contribution to the LCSR was computed in Ref. 01 



F?HQ 2 ) = 4 ^f«! CF (m\o) 2 , 



3/ 2 Q 4 



(3.8) 



by means of the quark condensate density. 



IV. CONSTRAINTS ON THE PION DAS 



The LCSR expression for the pion electromagnetic FF and the twist-4 DAs obtained in the framework of the 
renormalon approach can be used to extract constraints on the input parameters 62(^0) an< ^ ^4(^0)- ^ n order to 
perform numerical computations, we need to fix values of various parameters appearing in the relevant expressions. 
Namely , we take the Borel parameter M 2 within the interval 0.8 < M 2 < 1.5 GeV 2 and accept for the factorization 
and renormalization scales the following value: 



/4 



uQ 2 + uM 2 . 



(4.1) 



For the QCD coupling Q;s(Mr) the two-loop expression with A3 = 0.34 GeV is used. The value of the duality parameter 
so = 0.7 GeV 2 is borrowed from Shifman-Vainshtein-Zakharov sum rule |23q for the correlator of two u^^^d currents. 
The normalization scale is set equal to //§ = 1 GeV 2 . 

The Borel parameter dependence of the LCSR for different values of Q 2 is depicted in Fig. [21 From this figure one 
can conclude that the prediction for the FF is rather stable in the exploring range of M 2 . In what follows we choose 
the Borel parameter equal to M 2 = 1 GeV 2 . 

The scaled FF Q 2 F 7r (Q 2 ) and its different components are depicted in Fig. [21 In the calculations, the pion asymptotic 
DA and twist-4 DA ip£ (it) obtained from the renormalon approach are used. As is seen at Q 2 ~ 7.5 GeV 2 , the twist- 
4 contribution to the form factor exceeds the twist-2 one. This is important consequence of the higher conformal 
spin (renormalon) effects containing in the DA tpy' (u) . In Ref. |l9j the twist-4 term was calculated employing the 
asymptotic, i.e. the lowest conformal spin, form of g2Tr(u, This form leads to the combination 



20 

ip 4 (u, n 2 F ) = —5 2 {^ 2 F )uu [1 - u(7 - 8u)} , 



(4.2) 



with 



* 2 04) = * 2 (mu) 



«s(Mf) 



8Cf/3/3o 
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Q 2 (GeV 2 ) 



FIG. 3: The pion electromagnetic FF as a function of Q . The results are obtained employing the asymptotic DA. The solid 
line corresponds to the sum of the all contributions 13.51 . The dotted line shows O(qs) correction to the twist-2 term. 

^ 0,15- 




10 

Q 2 (GeV 2 ) 



FIG. 4: The twist-4 term as a function of Q . The curve 1 is computed using the ordinary asymptotic twist-4 DA. The 
twist-4 DAs obtained employing the renormalon method lead to the predictions shown by the line 2 and the broken lines. The 
correspondence between the lines and the parameters 62(^0), ^(Mo) ls: the line 2, 62(^0) = 0, ^(Mo) = 0; t ne dashed line, 
Mmo) = °- 17 > &4 0-*o) = -0.05; the dot-dashed line, 6 2 (mo) = 0.2, 64(^0) = 0. 



Comparing the twist-4 contributions found in the context of the different methods, one reveals their interesting 
features. The corresponding predictions are plotted in Fig. Here the curves 1 and 2 are computed using the 
standard asymptotic DA (see Eq. I|4.2|) 'l and the (p^\u) from the renormalon approach with 62(^0) = ^(Mo) = 0j 
respectively. The main difference between them is that the higher conformal spin effects shift the maximum of the 
twist-4 term towards larger values of Q 2 . This feature of the twist-4 term is more pronounced for DAs with &2(Mo) 7^ 0j 
^4(a*o) 7^ (the broken lines in Fig. 0J. Indeed, if the curve 2 takes its maximal value at Qq ~ 6 GeV 2 , for the broken 
lines we find Q\ ~ 11 GeV 2 . Starting from Qq, the twist-4 term slowly decreases, remaining larger than the standard 
prediction (the curve 1). Such a modification of the asymptotic behavior is another effect of the higher conformal 
spins. 

Calculations of Ref. correspond essentially to the "minimal" model of the twist-4 effects, where the restriction 

to the lowest conformal spin (a few lowest spins) probably underestimates the effect, while the renormalon model is 
a "maximal" model, where these effects are probably somewhat overestimated. Therefore, the renormalon model for 
the twist-4 DAs allows us, for the first time, to put a theoretically justified bound on the twist-4 contribution to the 
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FIG. 5: The la region for the scaled pion FF Q 2 F n (Q 2 ). The data arc taken from Refs. Q (the circles) and (the 
rectangles). In the la analysis only the solid data points are used. For the central solid curve the input parameters are 
fc^) =0.2, 6 4 (mo) = 0. 



pion form factor. Actually the change in absolute value of the twist-4 correction is not too dramatic, as one might 
expect. Thus, the ratio 



twA ren (Q 2 ) 
twA stand (Q 2 ) 

for the values b 2 {pl) = h(f4) = is c 1 ual to L25 at Q 2 = 10 GeV 2 and to 2 - 45 at Q 2 = 100 GeV 2 . 

In the renormalon approach, the twist-4 DAs and, hence, the twist-4 contribution to the form factor, depends on 
the input parameters and is not a constant background for the lea ding twist contribution. Therefore, performed la 
analysis results in conclusions, which differ from those made in Ref. [19|. In Fig. [3] we demonstrate the results of such 
la analysis. The data points included into the fitting procedure are shown by the solid points. Here we take into 
account the data Q 2 > 1.18 GeV 2 reported in Ref. [24|, and two new data points at Q 2 — 1; 1.6 GeV 2 obtained by 
the Fa collaboration |2^. From this analysis we extract the value of the input parameter 62(1 GeV 2 ), in the case of 
the DA with one nonasymptotic term 

b 2 (l GeV 2 ) = 0.2 ± 0.03. (4.3) 

For the pion DA with two nonasymptotic terms, we get 

b 2 (l GeV 2 ) = 0.2 ± 0.03, b 4 (l GeV 2 ) = -0.03 ± 0.06. (4.4) 

From Eqs. I|4.3|l and (|4.4|) it becomes evident that impact of ^2(^0) 011 the numerical computations is more important 
than a role played by 64(^0). Our analysis does not exclude also DAs with two positive input parameters. But 
it is worth noting that in our consideration we have used the data pij which were extracted indirectly from the 
pion electroproduction experiments through a model-dependent extrapolation to the pion pole. Moreover, the points 
Q 2 > 2 GeV 2 are imprecise suffering from the large errors and they are rather sparse. Correct and direct measurements 
of the form factor at Q 2 > 2 GeV 2 will improve the performed analysis and allow one to put more strong constraints 
on the pion DAs. 

In the region 1 GeV 2 < Q 2 < 10 GeV 2 , our prediction for the pion electromagnetic FF can be fitted to the following 
formula: 



(4.5) 
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FIG. 6: The pion leading twist DA extracted in this work. The scale Q 2 is fixed at 1 GeV 2 . For the central solid curve 
&2O0) = 0.2. For comparison the asymptotic distribution amplitude is also shown (dashed curve). 




FIG. 7: The two-particle twist-4 DAs ip^ (u, Q 2 ) (a), and — <^4 (u, Q 2 ) (b) obtained within the renormalon approach and 
using the constraint on the parameter b2(fx 2 ) 14.31 1. By dashed lines, for comparison, we plot the DAs obtained also within the 
renormalon approach, but using the pion asymptotic leading twist DA. 

where the uncertainties in the numerical coefficients are connected with the experimental errors. 

The pion twist-2 and two-particle twist-4 DAs calculated using the parameter l|4.3|l are shown in Figs. and [7| 
The shaded areas in the figures are obtained varying the parameter 62(1 GeV 2 ) within the allowed interval. The pion 
leading twist DA in the middle point uq = 0.5 takes the values 

ip {2) (0.5,1 GeV 2 ) = 1.05 T 0.07, tp {2) (0.5, 10 GeV 2 ) = 1.18 T 0.05. (4.6) 
This estimate is rather precise and does not contradict to the old Braun-Filyanov result, 

</? (2) (0.5) = 1.2 ±0.2, 

from the second paper in Ref. . The model DAs corresponding to Eq. (|4.4|l have the similar behavior and are not 
shown in the figures. 
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V. CONCLUDING REMARKS 

In the present work we have used the renormalon approach to determine the pion twist-4 DAs. In this approach 
the higher twist DAs are expressed in terms of the leading twist DA unambiguously. This fact has allowed us to 
avoid expansion of the higher twist DAs over the conformal spin and, at the same time, to take into account higher 
conformal spin effects. Of course, the renormalon approach is not suitable to model higher spin effects in the leading 
twist DA. Nevertheless, it considerably restricts a possible form of the higher twist DAs. 

The obtained model DAs have been employed for computation of the pion electromagnetic FF within the QCD 
LCSR method. For this purpose the correct expression for the twist-4 contribution to F 7r {Q 2 ) has been used [l9j 
and from comparison with the available experimental data the constraints on the input parameters 62(^0), ^4(^0) a ^ 
fj,Q = 1 GeV 2 have been deduced. 

The pion twist-2 DA <pir(u, Q 2 ) was an object of numerous investi gati ons. It was modeled using the various 
theoretical schemes and exclusive processes (see, for example, Refs. P. l2Ct l2^| ) . The models found in the present 
work are close to ones predicted in Ref. |20j . In Ref. [20( the power-suppressed corrections to F v (Q 2 ) were evaluated 
in the framework of the standard hard-scattering approximation and the running coupling method, which resulted 
in the Borel resummed FF [Q 2 F 7r (Q 2 )] res , whereas in the present paper we have computed the twist-4 contribution 
to F 7T (Q 2 ) in the context of the LCSR method and the renormalon-inspired twist-4 DAs. The new contribution of 
this work is that the renormalon approach has allowed to put an upper bound on the twist-4 contribution to the sum 
rules and obtain estimates, for the first time, of the effects due to higher conformal spins. We have gotten similar 
values of the pion DA parameters compared to other studies, so the LCSR approach seems to be protected from large 
uncertainties coming from higher twist corrections. 
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